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Instream  Flow  Studies 
Pennsylvania  and  Maryland 

Summary  Report 


EXECUTIVE  SUMMARY 

Existing  procedures  for  determining  mstream  flow  protection  levels  have  certain  deficiencies, 
which  result  in  conflicts  betw’een  agencies  that  regulate  water  supply  withdrawals  and  agencies  that  manage 
fisheries.  To  overcome  these  deficiencies,  the  Pennsylvama  Department  of  Environmental  Protection,  the 
Susquehanna  River  Basin  Commission,  Pennsylvania  Fish  and  Boat  Commission.  U.S.  Army  Corps  of 
Engineers,  Maryland  Department  of  the  Environment,  and  the  Biological  Resources  Division  of  the  U.S. 
Geological  Survey  cooperatively  conducted  an  instream  flow  needs  assessment  study.  The  goal  of  the 
study  is  to  develop  a procedure  for  determining  instream  flow  protection  levels  that:  (1)  is  based  on  tisher\' 
resouree  protection;  (2)  is  clearly  applicable  to  Pennsylvania  streams;  (3)  does  not  require  expensive  site- 
specific  studies;  and  (4)  can  be  easily  applied  during  the  administrative  review^  of  applications  for  surface 
water  allocations. 

The  basic  approach  to  the  development  of  the  procedure  is  to  conduct  mstream  flow  needs 
assessments  at  sites  selected  to  be  representative  of  a study  region,  and  then  regionalize  the  results  of  the 
site-speeific  assessments  to  develop  the  procedure.  Only  sites  w'ith  reproducing  trout  populations  and 
drainage  area  less  than  100  square  miles  were  included  m the  study. 

Physical  habitat  components  of  the  Instream  Flow  Ineremental  Methodology'  were  applied  to 
seleeted  study  sites  m the  Ridge  and  Valley  Freestone.  Ridge  and  Valley  Limestone,  Unglaciated  Plateaus, 
and  Piedmont  Upland  study  regions  in  Pennsylvania  and  Maryland.  The  evaluation  species  are  brook  and 
brown  trout.  Habitat  suitability  criteria  were  seleeted  from  the  literature,  and  tested  to  see  if  they 
adequately  represented  habitat  usage  on  Pennsylvania  streams.  These  critena  were  found  not  to  be 
applicable  to  Pennsylvania.  New  criteria  were  developed  from  the  data  collected  for  the  transferability 
study. 


Study  streams  were  selected  from  available  information,  and  divided  into  segments  based  on  length 
of  the  stream.  Study  sites  were  selected  near  the  midpomt  of  each  segment.  All  study  sites  had  good 
access,  reproducing  trout  populations,  and  good  water  quality.  Field  data  and  hydraulie  modeling  pro\  ided 
estimates  of  the  amount  of  habitat  available  within  a specified  range  of  flows.  The  amount  of  habitat 
available  for  all  life  stages  present  in  a defined  season  of  the  year  was  determined  for  that  range  of  flow  s. 

A computer  program  was  developed  to  estimate  the  effects  of  withdrawals  and  passb\'  flow  s on 
physical  microhabitat  and  availability  of  flow  for  withdrawals.  The  program  estimates  a number  of 
statistics  of  the  impact  for  vanous  combinations  of  withdrawal  and  passby  flow  for  any  project  site  in  the 
study  regions,  including  the  long-term  (average  annual)  impact.  This  computer  program  was  am  with 
many  combinations  of  species,  withdrawal  and  passby  flow  for  selected  stud\’  sites  within  a giv  en  class  of 
study  sites  (smdy  region,  segment  class)  to  estimate  the  average  annual  reduction  in  habitat  resulting  from 
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each  combination.  These  results  were  used  to  prepare  graphs  of  constant  habitat  impact,  and  the  percent  of 
time  that  water  supply  is  unavailable,  for  different  levels  of  impact. 

These  impact  curves  can  be  used  to  develop  statewide  policies  regarding  wiiich  impact  curve(s) 
should  be  used  to  establish  passby  flows.  They  also  can  be  used  to  determine  impact  of  a proposed 
withdrawal  at  any  site  in  these  study  regions.  These  curves  also  can  be  used  by  water  purveyors  to  analyze 
stream  intake  alternatives  that  meet  state  fishery  protection  levels  on  cold  water  streams  having  drainage 
areas  less  than  100  square  miles.  The  determination  of  which  impact  curve(s)  to  use  will  have  to  consider 
costs  both  to  the  environment  and  to  withdrawal  users.  Obviously,  the  curve  with  the  lowest  habitat  impact 
provides  the  greatest  protection  to  the  fishery’  habitat.  However,  as  the  degree  of  habitat  protection 
increases,  so  does  the  percent  of  time  that  withdrawals  cannot  be  made  because  of  flow  limitations  or 
passby  flow  requirements. 

Although  regional  critena  have  been  developed,  the  computer  program  also  can  be  used  to  evaluate 
conditions  not  considered  in  the  development  of  the  regional  critena.  A regional  hydrology  procedure  has 
been  developed  to  provide  hydrology  for  the  computer  program. 

A detailed  description  of  the  methodology’  developed  and  applied  in  this  study,  and 
recommendations  for  additional  studies,  are  presented. 


1.0  INTRODUCTION 

This  report  is  a summary’  of  Instream  Flow  Studies — Pennsylvama  and  Maryland  (Denslinger  and 
others,  1998).  It  is  intended  to  summarize  the  procedures  and  the  results  of  the  instream  flow  needs 
assessment  study  conducted  in  Pennsylvania  and  Maryland.  The  purpose  of  the  study  is  to  develop  a 
regional  procedure  for  estimating  the  effects  of  withdrawals  and  passby  flows  on  fishery  habitat,  and  to 
develop  regional  cnteria  for  determining  passby  flows  for  use  during  the  administrative  review  of 
applications  for  water  supply  allocations. 

The  complete  report  is  available  from  the  Susquehanna  River  Basin  Commission  (SRBC). 


2.0  STUDY  SUMMARY 


2.1  Study  Need 

Historically,  the  Pennsylvania  Department  of  Environmental  Protection  (Pa.  DEP)  has  used  several 
different  procedures  to  determine  mandated  conservation  releases  from  major  water  supply  reservoirs,  or 
mandated  passby  flows  at  smaller  dams  and  intake  structures,  to  reduce  impacts  on  fishery  resources. 
These  procedures  do  not  directly  consider  the  effects  of  withdrawals  on  fishery  resources,  and  the 
conservation  flows  are  derived  from  hydrologic  records  utilizing  a statistical  low  flow.  In  1992,  Pa.  DEP 
and  the  Pennsylvania  Fish  and  Boat  Commission  (PFBC)  agreed  to  use,  to  the  extent  practicable,  a 
procedure  based  upon  the  Tennant  method  (Tennant,  1976).  The  method  attempts  to  incorporate  aquatic 
resource  needs,  but  has  questionable  application  to  Pennsylvama  streams,  and  may  unnecessarily  reduce 
the  yield  that  can  be  obtained  from  water  supply  sources,  while  providing  more  than  adequate  protection  to 
aquatic  resources. 
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To  correct  the  deficiencies  of  existing  methods.  Pa.  DEP,  Susquehanna  River  Basin  Commission 
(SRBC),  PFBC.  U S.  Army  Corps  of  Engineers  (US  COE).  Mary  land  Department  of  Environment  (MDE). 
Chesapeake  Bay  Program  (CBP),  and  U S.  Geological  Survey  Biological  Resources  Division  (GS3RD) 
conducted  this  study  to  develop  a procedure  for  determining  mstream  flow  protection  levels  that:  (1)  is 
based  on  fisher\’  resource  protection:  (2)  is  clearly  applicable  to  Pennsylvania  streams:  (3)  does  not  require 
expensive  site-specific  studies;  and  (4)  can  be  easily  applied  dunng  the  administrative  review  of 
applications  for  surface  water  allocations. 

SRBC’s  responsibilities  for  managmg  the  water  resources  of  the  Susquehanna  basin  include 
protecting  mstream  flows  tlirough  the  regulation  of  ( 1 ) certain  water  withdrawals  where  signatories  to  the 
Susquehanna  Compact  (Susquehanna  River  Basin  Commission,  1972)  do  not  have  the  authority:  and  (2) 
consumptive  use  of  water. 

The  State  of  Maiyland  also  is  interested  in  the  development  of  new  procedures  for  determining 
flow's  that  protect  the  biota  and  also  allow  water  supply  withdrawals.  The  Maiyland  water  suppK 
allocation  program  uses  the  Maryland  Most  Common  Flow  Method  to  establish  conservation  flow 
requirements  for  water  supply  withdrawals  and  reservoir  projects. 

There  are  many  important  instream  flow  protection  issues.  Among  the  pnoritv'  issues  are; 

• The  effects  of  withdrawals  and  consumptive  uses  on  aquatic  biota  in  cold  water  trout  streams: 

• The  effects  of  withdrawals  and  consumptive  uses  on  aquatic  biota  in  tnbutaiy  streams  with 
warm  w'ater  fisheries; 

• The  effect  of  withdrawals  and  consumptive  uses  on  the  aquatic  biota  in  major  rivers:  and 

• The  effect  of  consumptive  uses  on  the  receiving  waters  of  the  Chesapeake  Bay. 

The  issue  of  the  effect  of  withdrawals  and  consumptive  uses  on  aquatic  biota  in  cold  water  streams 
is  considered  the  most  important  because  of  existmg  critical  conflicts  between  withdrawals  and  mstream 
uses  on  those  streams. 

2.2  Study  Purpose  and  Approach 

The  purpose  of  this  study  is  to  develop  a procedure  for  determining  instream  flow  needs  for 
streams  with  naturally  reproducing  trout  populations  (reproducing  trout  streams)  in  portions  of 
Pennsylvania  and  Maryland,  that  does  not  require  a stream-specific  impact  analysis  study.  The  procedure 
must  be  based  on  fishery  habitat,  and  instream  flow  needs  must  be  derived  from  hydrologic  data  and  the 
data  developed  in  the  study. 

The  basic  approach  to  the  problem  is  to  conduct  instream  flow  needs  assessment  studies  at  selected 
representative  sites,  and  then  regionalize  the  results  of  the  site-specific  assessments  to  develop  the 
procedure.  Because  of  existing  cntical  conflicts  between  instream  and  withdrawal  uses  on  small  trout 
streams  in  the  Ridge  and  Valley,  and  unglaciated  parts  of  the  Appalachian  Plateaus.  ph\'siographic 
provinces,  this  study  focuses  on  those  areas.  Some  streams  in  the  Piedmont  Upland  physiographic  section 
in  Maiy'land  also  were  studied.  Only  reproducing  trout  streams  with  drainage  areas  less  than  100  square 
miles  are  included  in  this  study,  because  most  reproducing  trout  streams  in  Pennsylvania  and  Marv  land  are 
in  that  size  range,  and  because  the  effects  of  withdrawals  on  instream  uses  are  most  critical  on  those 
streams.  Tlierefore.  the  procedure  applies  onh'  to  those  streams  at  this  time. 


2.3  Instream  Flow  Needs  Assessment  Methodologies 

The  Instream  Flow  Incremental  Methodology  (IFIM)  (Bovee,  1982)  and  the  wetted  perimeter 
method  (Collings.  1974;  Nelson,  1984;  Leathe  and  Nelson,  1989)  were  both  applied  to  selected  streams  in 
this  study.  IFIM  is  the  most  sophisticated  method  available  for  determinmg  instream  flow  needs,  and  is 
specifically  designed  to  assess  effects  of  man-made  changes  m flow  on  the  habitat  available  for  fish.  The 
wetted  perimeter  method  has  been  used  by  other  mvestigators  to  establish  instream  flow  protection  levels. 

The  ifiM  methodology  includes  physical  microhabitat  components  shown  in  Figure  2.1,  which 
were  used  in  this  study  to  estimate  impacts  of  different  combinations  of  natural  flows  and  withdrawals. 
The  Physical  Habitat  Simulation  (PHABSIM)  computer  programs  were  used  to  evaluate  physical 
microhabitat.  The  methodology’  uses  evaluation  species  selected  for  recreational,  economic,  or  ecological 
importance  to  evaluate  impacts  of  withdrawals  on  stream  ecology.  Habitat  suitability  criteria  (HSC)  are 
used  to  describe  the  usability  of  depth,  velocity,  substrate,  and  cover  for  each  life  history  stage  (adult, 
juvenile,  fry  and  spawning)  for  each  evaluation  species.  Depth  measurements  are  made  at  different  flows 
for  a number  of  measurement  points  across  one  or  more  transects  at  a study  site.  Velocity,  substrate  and 
cover  measurements  are  made  at  the  same  measurement  points  at  one  flow.  The  depth  and  velocity 
measurements  are  used  to  calibrate  hydraulic  models  for  each  measurement  point.  The  calibrated  models 
are  used  to  simulate  depth  and  velocity  at  a number  of  flows  within  a specified  range,  for  each  transect  and 
measurement  point.  The  simulated  depth  and  velocity  values  are  combined  with  substrate  and  cover 
measurements,  and  with  the  HSC,  to  generate  a relationship  between  habitat  and  flow  for  each  transect, 
species  and  life  stage.  Habitat  is  defined  as  weighted  usable  area  (WUA).  The  difference  between  WUA 
available  at  different  flows  with  and  without  a proposed  project  is  the  impact  of  the  project  on  the  habitat 
available. 

The  wetted  perimeter  of  a transect  generally  increases  rapidly  with  flow  for  flows  less  than  some 
amount,  and  increases  less  rapidly  for  higher  flows.  Wetted  perimeter  is  plotted  versus  flow,  and  the  flow 
value  at  the  point  where  the  graph  changes  slope  (inflection  point)  is  assumed  to  be  the  amount  of  flow 
necessary  to  protect  the  biota.  The  wetted  perimeter  method  is  usually  applied  to  the  riffle  area(s)  of 
streams,  because  those  are  the  most  productive  areas  for  aquatic  invertebrates,  which  are  the  food  base  for 
certain  species  of  fish  (Collings,  1974;  Nelson,  1984;  Leathe  and  Nelson,  1989). 

The  wetted  perimeter  method  has  the  advantage  of  being  quick  and  inexpensive  to  apply. 
However,  the  method  has  a number  of  questionable  assumptions  and  limitations  (Leonard  and  others,  1986; 
Mohrhardt,  1987).  The  major  assumption  is  that  the  flow  at  the  inflection  point  needs  to  be  maintained  to 
ensure  an  adequate  food  supply  for  the  fish,  but  this  assumption  has  not  been  verified.  The  method  does 
not  allow  evaluation  of  the  effects  of  withdrawals  on  the  biota. 

Leathe  and  Nelson  (1989)  list  five  major  factors  (Hall  and  Knight,  1981)  that  control  fish 
abundance  in  streams:  streamflow;  habitat  quality;  food  abundance;  predation;  and  movement  and 
migration.  Any  of  these  may  be  the  limiting  factor  for  any  given  stream.  Standard  setting  methods  such  as 
the  wetted  perimeter  method  identify  minimum  flow  standards,  while  incremental  methods  such  as  IFIM 
quantify  tradeoffs  between  withdrawals  and  instream  uses  by  examining  the  response  of  fish  habitat  to 
changes  in  flow  (Leathe  and  Nelson,  1989). 

Wliere  man-made  changes  in  streamflow  such  as  withdrawals  limit  the  amount  of  habitat  available, 
a method  that  evaluates  the  effects  of  incremental  changes  in  streamflow  such  as  IFIM,  is  probably  most 
appropriate.  The  wetted  perimeter  method  may  be  appropriate  where  food  supply  is  the  limitmg  factor,  or 
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Real  World  River  System 


Figure  2. 1.  Components  of  Instream  Flow  Incremental  Methodology 


when  a simple  method  is  needed  to  develop  basinwide  standards  for  use  in  preliminary  watershed  planning 
(Leathe  and  Nelson.  1989). 

2.4  Evaluation  Species  and  Habitat  Suitability  Criteria 

Brook  and  brown  trout  were  selected  as  representative  species  for  the  evaluation  of  habitat 
availability  and  the  impact  of  withdrawals.  These  species  were  selected  because  they  are  the  most 
important  economically  and  recreationally  in  the  study  regions.  The  periods  when  the  different  life  stages 
of  these  species  are  present  were  determined  and  used  to  define  seasons  for  impact  analysis. 

For  these  species,  depth  and  velocity'  suitability'  criteria  were  selected  from  the  literature 
(Bovee,  1978;  Bovee,  oral  communication,  1994;  Aceituno  and  others,  1985;  Raleigh  and  others,  1986; 
Jirka  and  Homa,  1990;  Harris  and  others,  1992;  Normandeau  Associates  Inc.,  1992;  and  Gary'  Whelan. 
Michigan  Department  of  Natural  Resources,  oral  commumcation,  1994)  for  use  in  the  IFIM  methodology. 
The  Normandeau  (1992)  criteria  for  adult  and  juvenile  life  stages,  the  Bovee  (1978)  criteria  for  the  fry  life 
stage,  and  the  Whelan  (oral  communication,  1994)  criteria  for  the  spawning  life  stage  were  selected  for 
testing.  A substrate/cover  classification  scheme  and  correspondmg  suitability  criteria  were  developed  for 
use  in  the  study,  based  on  professional  judgement. 

In  accordance  with  GSBRD  recommendations,  these  suitability  criteria  were  tested  to  determine 
whether  they  could  be  transferred  to  Permsylvania.  Four  streams  were  selected  for  transferability  testing, 
one  brook  trout  stream,  and  one  brown  trout  stream,  in  both  the  Ridge  and  Valley  Freestone  and  the 
Unglaciated  Plateau  study  regions,  respectively.  The  transferability  study  generally  followed  the 
methodology  described  by  Thomas  and  Bovee  (1993).  Depth,  velocity,  substrate  and  cover  were  reeorded 
at  locations  occupied  by  the  various  life  stages  of  the  evaluation  species,  and  at  locations  not  occupied  by 
fish.  Statistical  analyses  of  these  data  showed  that  the  selected  criteria  were  not  suitable  for  use  in 
Pennsylvania.  New  suitability  criteria  were  developed  using  the  data  collected  for  the  transferability  study, 
and  used  in  the  subsequent  PHABSIM  studies. 

2.5  Study  Regions  and  Study  Stream  Selection 

To  develop  a regional  method,  reproducing  trout  streams  were  classified  according  to  key  physical 
features  that  have  a direct  influence  on  the  physical  variables  and  stream  attributes  used  to  quantify  fishery 
habitat.  Streams  were  classified  according  to  study  region,  species,  and  segment  number. 

Study  regions  were  based  on  physiographic  provinces  and  sections  (Fenneman,  1938;  Pennsylvania 
Dept,  of  Environmental  Resources,  1989).  In  the  Ridge  and  Valley  physiographic  province,  streams  were 
classified  into  study  regions  based  on  limestone  (including  dolomite)  or  freestone  (e.g.,  sandstone,  shale, 
conglomerate)  geology,  rather  than  physiographic  sections.  In  the  Appalachian  Plateaus  physiographic 
province,  streams  were  classified  into  glaciated  and  unglaciated  study  regions,  based  on  the  location  of  the 
glacial  boundary  (Pennsylvania  Dept,  of  Environmental  Resources,  1989).  Streams  in  the  unglaciated 
physiographic  sections  were  combmed  into  one  study  region  called  the  Unglaciated  Plateau.  Trout  streams 
in  the  Piedmont  Province  were  classified  based  on  physiographic  section  and  limestone/freestone  geology. 
Physiographic  regions  and  sections,  and  parts  of  the  Ridge  and  Valley  Province  underlain  by  limestone,  are 
shown  in  Plate  1 . 

In  Pennsylvania,  the  conflicts  between  withdrawal  and  instream  uses  are  most  critical  on  cold 
water  streams  in  the  Ridge  and  Valley  physiographic  province,  and  the  unglaciated  part  of  the  Appalachian 
Plateaus  Province.  The  Piedmont  Upland  physiographic  section  was  selected  because  there  are  more  trout 
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streams  m that  section  than  other  sections  in  Mary  land.  For  these  reasons,  and  because  of  time  and  cost 
constraints,  only  the  Ridge  and  Valley  Limestone.  Ridge  and  Valley  Freestone.  Unglaciated  Plateau,  and 
Piedmont  Upland  (freestone)  study  regions  were  included  in  this  sTud\'. 

Parts  of  five  counties  in  the  southwestern  comer  of  Pennsylvania  were  deleted  from  the  stud\ 
because  the  streams  have  very  low  yield,  and  there  are  few  reproducing  trout  streams  m the  area.  For  these 
reasons,  there  are  few  water  supply  withdrawals  from  reproducing  trout  streams  in  that  area,  Tlie  counties 
deleted  are  Beaver,  Allegheny,  Washington,  Greene,  and  part  of  Fayette. 

Lists  of  streams  with  naturally  reproducing  trout  populations  (reproducing  trout  streams)  in  each 
studv  region  in  Pennsylvania  were  developed  from  existing  PFBC  and  Pa.  DEP  data.  The  list  of  trout 
streams  in  Maryland  was  developed  from  a report  prepared  by  Steinfelt  (1991).  The  presence  of 
reproducing  trout  populations  on  certam  study  streams  selected  m Pennsy  K ama  was  venfied  in  the  field 
because  the  PFBC  records  were  incomplete.  Potential  study  streams  were  selected  from  these  lists  by 
stratified  random  sampling.  The  actual  study  streams  were  selected  m the  field  from  the  list  of  potential 
streams,  also  using  stratified  random  sampling.  All  study  streams  had  good  access,  reproducing  trout 
populations,  good  water  quality,  and  no  significant  human  influences. 

Study  streams  were  divided  into  segments  based  on  stream  length,  which  was  used  as  a surrogate 
for  stream  slope.  The  maximum  allowable  length  of  stream  segments  was  set  at  five  miles,  based  on 
statistical  analysis  of  stream  length  data.  The  actual  segment  length  depended  on  the  total  length  of  stream. 

The  final  selected  study  sites  are  shown  in  Plate  1. 

A key  assumption  is  that  a total  of  30  study  segments  is  adequate  to  represent  the  variability'  in 
hydrology  and  habitat  response  to  withdrawals  m each  study  region.  Approximately  30  segments  of 
various  size  were  studied  in  each  of  the  three  study  regions  in  Pennsylvania  (Ridge  and  Valley  Limestone. 
Ridge  and  Valley  Freestone,  Unglaciated  Plateau),  but  due  to  time  and  cost  constraints,  only  12  segments 
were  studied  in  the  Piedmont  Upland  study  region.  The  proportion  of  streams  in  each  segment  class  w as 
approximately  equal  to  the  proportion  of  streams  in  that  class  in  the  entire  number  of  reproducing  trout 
streams  in  the  respective  region. 

2.6  Field  Data  Collection 

Once  the  study  streams  were  chosen,  a study  site  was  selected  near  the  midpoint  of  each  segment. 
Then  the  relative  amount  of  each  different  mesohabitat  type  (nffle.  run.  pool)  was  estimated  for  each  study 
site.  A representative  occurrence  of  each  mesohabitat  type  was  selected,  and  transects  were  located  near 
the  midpoint  of  the  respective  mesohabitat  type. 

Flow  rate  and  water  surface  elevation  were  measured  at  each  transect,  at  a sufficient  number  of 
flows  to  allow  calibration  of  a hydraulic  model.  The  range  of  measurement  flows  was  specified  to  satisfv 
extrapolation  enteria  for  hydraulic  calibration  (U.S.  Geological  Survey,  Biological  Resources  Di\ision. 
1994),  and  simulate  flows  over  the  range  beKveen  maximum  and  mimmum  median  monthly  flow  s at  each 
site.  Velocity  distnbution,  substrate  and  cover  were  measured  at  a number  of  points  across  each  transect, 
generally  at  only  one  flow.  The  measurement  points  were  selected  to  represent  changes  in  habitat  or 
velocity^  across  the  transect.  Field  data  collection  followed  standard  procedures  (Bovee.  undated: 
Buchanan  and  Somers.  1969). 
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Data  were  collected  to  show  the  location  of  trout  redds  (nests)  within  each  mesohabitat  type,  to 
evaluate  whether  transects  located  near  the  midpoint  of  a mesohabitat  type  adequately  represented 
spawning  habitat.  For  each  mesohabitat  tspe,  a large  proportion  of  the  redds  were  found  in  the  central  half 
of  that  tvpe.  Tlierefore,  it  was  concluded  that  transects  located  near  the  midpoint  of  each  mesohabitat  type 
adequately  represented  spawning  habitat. 

The  following  problems  were  encountered  during  the  study  site  selection  and  field  data  collection 
phases  of  the  study.  In  limestone  regions,  aquatic  vegetation  frequently  caused  difficulty  in  obtaining  valid 
velocit\'  and  flow  measurements.  Seasonal  changes  in  vegetation  resulted  in  changes  in  depth,  velocity  and 
roughness  for  different  measurements,  which  made  hydraulic  calibration  difficult,  and  in  some  cases 
impossible.  For  some  streams,  changes  in  transect  geometry  between  measurements,  usually  as  a result  of 
high  flows,  also  caused  inconsistencies  between  measurements  and  required  collection  of  additional  data,  or 
m some  cases,  deletion  of  a study  site.  Streams  in  the  Piedmont  Upland  in  Maryland  showed  signs  of 
unstable  bed  and  banks.  Future  hydraulic  and  habitat  conditions  for  these  streams  may  be  different  from 
current  conditions,  so  the  habitat  analyses  should  be  used  with  caution. 

2.7  Hydrology  and  Habitat  Modeling 

Hydrology'  was  developed  from  flow  data  collected  at  stream  gages  selected  to  be  representative  of 
the  study  streams.  The  follovUng  hydrology  was  developed  for  the  study  sites: 

• Average  daily  flow  (ADF); 

• Median  flow  for  the  entire  period  of  record; 

• Median  monthly  flows  for  each  month  for  the  entire  period  of  record; 

• Time  senes  of  median  monthly  flows;  and 

• Annual  and  seasonal  flow  duration. 

The  hydrology  for  each  study  site  was  generally  developed  from  the  corresponding  hydrology  for  a 
selected  stream  gage,  by  multiplying  flows  at  the  gage  by  the  ratio  of  drainage  area  at  the  site  to  drainage 
area  at  the  gage.  Stream  gages  were  selected,  based  on  drainage  area  size,  proximity  to  the  study  site, 
similar  geology  and  topography,  and  judgment.  For  study  streams  or  gages  with  mixed  limestone  and 
freestone  geology,  significant  springs,  withdrawals,  or  wastewater  treatment  plant  flows,  more  complex 
procedures  were  used  to  derive  the  hydrology. 

Hydraulic  models,  based  on  Manning's  equation  (Chow,  1959;  Bovee,  1982),  were  calibrated  for 
each  transect  and  measurement  point  at  each  study  site.  The  calibrated  hydraulic  model  was  used  to 
simulate  velocity  and  depth  for  18  flows  in  the  range  between  maximum  and  minimum  median  monthly 
flows,  in  accordance  with  established  extrapolation  cnteria  (U.S.  Geological  Survey,  Biological  Resources 
Division.  1994).  The  simulated  depth  and  velocity  data  were  combined  with  the  substrate  and  cover  data, 
and  the  HSC,  to  develop  WUA  versus  flow  relationships  for  each  evaluation  species,  life  stage,  and 
transect.  Tlie  percentages  of  each  mesohabitat  type  for  each  study  site  were  used  to  compute  a weighted 
average  WUA  versus  flow  relationship  for  each  study  site,  species,  and  life  stage. 

A pilot  study  was  conducted  to  determine  whether  binary  HSC  should  be  used  instead  of  univariate 
HSC.  as  recommended  by  Bovee  and  others  (1994).  Univanate  criteria  can  have  values  over  the  entire 
range  from  1 to  0,  but  binary  cntena  can  only  have  a value  of  either  1 or  0.  In  effect,  binary  criteria 
discard  habitat  that  may  be  usable  by  the  fish  species,  but  that  is  not  optimal  for  those  species.  The 
univariate  cntena  developed  from  the  transferability’  study  data  were  modified  to  binary  form. 
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WL^A  versus  flow  relationships  for  each  t\pe  of  criteria  were  computed  and  plotted.  The 
univariate  cntena  resulted  in  smooth  WUA  plots  for  all  streams  considered.  The  binarv  cnteria  resulted  in 
highlv  vanable.  sometunes  saw-tooth  WUA  plots,  with  much  lower  amounts  of  habitat,  which  made 
interpretation  of  these  plots  difficult.  The  WUA  curves  based  on  univanate  cnteria  appeared  more  realistic 
and  consistent  with  expected  relationships  for  the  study  streams,  which  support  good  trout  populations. 
The  marginal  habitat,  which  is  not  considered  in  tlie  binary  cntena.  may  be  vciy  imponant  to  the  trout 
populations.  For  that  reason,  univanate  cntena  were  used  to  develop  the  WUA  relationships  used  in  the 
impact  analysis. 

2.8  Wetted  Perimeter  Analysis 

Wetted  penmeter  versus  flow  plots  were  prepared  using  the  output  from  the  hydraulic  simulations, 
for  the  riffie  transects  only.  This  procedure  effectively  assumes  the  inflection  point  (change  of  slope)  of  the 
plot  occurs  in  the  range  between  maximum  and  minimum  median  monthly  flow.  The  flow  rates  at  the 
inflection  points  of  the  curves  were  tabulated  for  each  study  region,  and  converted  to  flow  rates  per  unit 
area  and  to  percent  of  ADF.  These  tabulations  showed  a lot  of  variabilitv’  of  the  flow  rates  at  the  inflection 
points  within  each  region. 

The  plots  were  extrapolated  to  zero  wetted  perimeter  at  zero  flow.  The  extrapolation  substantially 
changed  many  graphs,  and  usually  introduced  a lower  inflection  point.  The  resulting  inflection  points  also 
were  tabulated  for  the  three  study  regions  in  Pennsylvania,  and  were  generally  lower  than  the  inflection 
points  determined  from  the  simulation  flows  alone.  The  conclusion  is  the  wetted  penmeter  data  developed 
from  the  limited  range  of  simulation  flows  are  not  adequate  to  allow  selection  of  inflection  points. 
Therefore,  comparisons  with  the  results  of  the  IFIM  method  are  not  possible  without  collecting  additional 
extreme  low  flow  data. 

2.9.  Impact  Assessment  Methods  and  Results 

The  median  monthly  habitat  was  assumed  to  be  the  best  measure  of  the  amount  of  habitat  t\picallv 
available.  A pilot  study  showed  that  the  habitat  available  at  the  median  monthly  flow  is  essentially  the 
same  as  the  median  of  the  daily  habitat  determined  from  daily  flows.  Therefore,  the  median  monthly 
habitat  was  defined  as  the  habitat  value  associated  with  the  median  monthly  flow  for  subsequent  anah  ses. 

To  obtain  WUA  versus  flow  relationships  for  each  study  site,  each  species,  and  each  season,  the 
life  stage  present  in  that  season  with  the  least  habitat  at  any  flow  was  assumed  to  be  the  most  critical  life 
stage  to  be  protected  at  that  flow  (Orth  and  Leonard.  1990).  A procedure  was  developed  and  unplemented 
to  compute  these  relationships,  which  are  called  renormalized  minimum  weighted  usable  area  (RiVfWUA). 

To  determme  a conservation  flow  that  would  protect  the  habitat  available,  two  alternative 
definitions  of  habitat  loss  were  considered,  no-loss  of  habitat,  and  no-net-loss  of  median  monthh'  habitat. 
For  this  study,  no-loss  of  habitat  was  defined  as  no  reduction  in  RJVfWUA  at  any  flow.  No-net-loss  of 
habitat  was  defined  as  no  reduction  of  RjVfVVUA  at  the  median  monthly  flow.  The  no-loss  cnterion 
unnecessanly  limits  the  withdrawals  under  a wide  range  of  conditions,  considering  that  natural  flow  and 
available  habitat  fluctuate  within  months,  and  \ ears.  and  among  years.  Analysis  of  the  no-net-loss  critenon 
for  1 1 study  sites,  for  both  brook  trout  and  brown  trout,  for  three  months  of  the  summer  season,  showed 
that  the  no-net-loss  flow  was  equal  to  the  median  monthly  flow  for  virtualK’  all  cases.  The  conclusion  was 
this  cnterion  severely  lunits  the  withdrawals  durmg  the  summer  season.  Therefore,  more  detailed 
procedures  were  developed  to  assess  the  impact  of  water  withdrawals  on  the  habitat  available. 


9 


The  purpose  of  impact  analysis  is  to  determme  the  magnitude  of  the  impact  of  withdrawals  on 
habitat  over  a full  range  of  flows,  and  to  use  that  information  to  establish  criteria  for  passby  flows.  The 
impact  IS  defined  as  the  percentage  difference  beuveen  habitat  (RJVIWUA)  available  without  the 
withdrawal,  and  the  habitat  available  with  the  withdrawal  in  place. 

Two  alternative  procedures  were  developed  to  estimate  the  impact  of  withdrawals  on  available 
habitat.  The  first  procedure  analyzes  the  effect  of  withdrawals  on  time  series  of  median  monthly  flow  and 
habitat.  The  second  procedure  analyzes  the  effect  of  withdrawals  on  flow  and  associated  habitat  duration. 

The  time  senes  impact  analysis  procedure  is  designed  to  estimate  the  long-term  effect  of 
withdrawals,  for  a specific  project  site  and  a specific  combination  of  withdrawal  and  passby  flow  using 
median  monthly  flow  time  series.  A time  series  is  simply  a set  of  values  arranged  in  chronological  order. 
The  method  also  can  be  used  with  other  time  steps  such  as  daily,  but  for  shorter  periods  of  record.  The 
procedure  estimates  the  average  regional  impact  at  a project  site,  in  a given  segment  class,  of  a combination 
of  withdrawal  and  passby  flow  (both  expressed  as  percentage  of  ADF)  by  determining  impacts  on  each 
study  site  in  that  class.  Then  the  impacts  are  averaged  across  the  study  sites  in  that  class. 

A computer  program  has  been  developed  in  Microsoft  Excel  7.0  format,  to  estimate  the  impact  of 
withdrawals  for  any  location  within  a study  region.  There  are  two  separate,  but  related,  computer 
programs  included  in  the  package.  The  first,  called  the  "detailed  analysis”  program,  provides  detailed 
estimates  of  the  average  impacts  of  any  combination  of  withdrawal  and  passby  flow  on  the  flow  and  habitat 
at  a project  site,  and  can  analyze  a number  of  different  combinations  of  species  and  trout  management 
procedures.  This  program  also  allows  determmation  of  the  percent  of  time  the  withdrawal  is  not  available 
for  the  given  combination  of  withdrawal  and  passby  flow.  The  second  program,  called  the  “preliminary 
analysis”  program,  provides  general  estimates  of  the  effect  of  withdrawals  and  passby  flows.  The  results 
of  the  prelimmaiy'  analysis  program  can  serve  as  a starting  point  for  more  complete  analysis  using  the 
detailed  analysis  program. 

The  detailed  analysis  program  has  been  used  with  the  hydrology  and  RMWUA  data  for  selected 
study  sites  to  develop  habitat  impact  curves  for  the  Unglaciated  Plateau,  Ridge  and  Valley  Freestone,  and 
Ridge  and  Valley  Limestone  study  regions.  One  study  site  was  selected  for  impact  analysis  to  represent 
each  stream  gage  used  to  develop  hydrology  for  each  segment  class  and  study  region.  The  data  for  the 
curves  were  obtained  by  systematically  varying  withdrawals  and  passby  flows.  For  each  segment  class  in 
each  region,  twenty-seven  combinations  of  withdrawal  and  passby  flows  (e.g.,  10  percent  ADF  withdrawal 
and  5 percent  ADF  passby  flow)  were  run  for  each  of  the  stream  gages  represented,  and  for  each  species 
considered.  Three  species  were  analyzed;  wild  brook  trout,  wild  brown  trout,  and  combined  wild  brook  and 
brown  trout. 

For  each  study  site,  the  average  annual  percent  reduction  in  RMWUA  across  the  period  of  record 
was  used  as  the  measure  of  impact.  Curves  of  constant  impact  (e.g.,  25  percent  reduction  of  habitat)  were 
developed  for  each  region,  species,  withdrawal,  and  passby  flow.  The  Ridge  and  Valley  Limestone  region 
was  split  into  two  groups,  based  on  whether  the  amount  of  limestone  on  the  watershed  was  greater  or  less 
than  50  percent,  and  different  curves  were  developed  for  each  group. 

Comparison  of  the  average  annual  impacts  for  the  selected  study  streams  within  each  region 
showed  little  variability  between  the  average  impacts  across  streams  and  the  maximum  and  minimum 
values  of  those  impacts.  This  companson  indicated  that,  while  hydrology  and  stream  characteristics  were 
highlv'  variable,  habitat  impacts  were  fairly  consistent  within  each  region.  However,  impacts  for  a given 
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combination  of  species,  withdrawal,  and  passby  flow  were  vers-  different  among  regions.  Tins  supported 
the  basic  study  concept  that  streams  would  react  similarly  within  regions,  but  differently  among  regions. 

For  the  Ridge  and  Valley  Freestone  study  region,  the  impact  curves  for  segment  classes  1.  2.  and  3 
were  close  together,  so  these  curves  were  averaged.  Because  segment  class  4 included  only  one  stream,  no 
impact  curves  were  provided  for  that  class. 

For  the  Ridge  and  Valley  Limestone  study  region,  the  average  annual  impacts  showed  significant 
scatter  among  streams.  These  study  sites  were  further  classified  based  on  the  percentage  of  limestone  in 
the  watershed,  which  significantly  reduced  the  scatter,  but  also  reduced  the  sample  size.  especialK'  for 
segment  classes  2,  3,  and  4.  Because  of  limited  sample  size  and  the  effect  of  existing  withdrawals.  WWTP 
flows,  and  springs  (or  caves)  on  the  hydrology'  at  these  study  sites,  impact  curves  were  developed  only  for 
segment  class  1 sites  in  this  study  region. 

A partial  list  of  limestone  streams  has  been  provided.  This  list  may  be  incomplete  and  should  be 
used  with  caution.  Streams  not  included  m the  list  should  be  classified  as  limestone  on  case-by-case  basis. 

In  the  Unglaciated  Plateau  study  region  Comparison  of  the  impact  curves  showed  a difference 
between  segment  class  1 and  class  2 sites.  There  were  no  segment  class  3 or  4 study  sites  in  that  region. 

The  final  constant  habitat  impact  curves  are  shown  in  Figures  2.2  through  2.11.  The  habitat 
impact  curves  for  Ridge  and  Valley  Limestone  Group  1 should  be  used  for  project  streams  and  sites  where 
more  than  50  percent  of  the  watershed  upstream  from  the  study  site  is  underlain  by  limestone,  and  the 
curves  for  Group  2 should  be  used  where  less  than  50  percent  of  the  watershed  is  underlain  by  limestone. 
For  all  three  study  regions,  the  impact  curves  for  brown  trout  and  combined  brown  and  brook  trout  are 
similar.  For  that  reason,  only  one  set  of  curves  is  shown  for  those  species  for  each  study  region.  There  are 
sigmficant  differences  between  impacts  on  brook  and  brown  trout,  as  well  as  between  brook  trout  and 
combined  brook  and  brown  trout. 

The  maximum  and  the  90  percent  probability'  of  exceedance  measures  of  habitat  impact  also  were 
considered.  The  average  impact  curves  show  the  long-term  effect,  and  maximum  impact  curves  show  the 
short-term  effect.  The  impact  curves  based  on  the  average  impact  are  included  in  this  report,  based  on  the 
assumption  that  long-term  average  impacts  to  habitat  may  result  in  average  impacts  to  fish  biomass  of 
similar  magnitude.  However,  since  short-term  maximum  impacts  to  habitat  may  have  more  acute  effects, 
both  long-term  and  short-term  impacts  should  be  considered  when  makmg  decisions  regarding  habitat 
protection.  The  impacts  at  the  90  percent  probability'  of  exceedance  were  found  to  be  very  close  to  the 
maximum  impacts,  and  thus  provided  no  advantage. 

The  constant-habitat-impact  graphs  (Figures  2.2  through  2.11)  also  show  the  impact  of  a given 
passby  flow  on  the  percentage  of  time  that  a given  withdrawal  is  not  available.  Obviously,  the  curve  with 
the  lowest  habitat  impact  provides  the  greatest  protection  to  the  fishery  habitat.  However,  as  the  degree  of 
protection  increases,  so  does  the  percent  of  time  that  withdrawals  cannot  be  made  because  of  passbv 
requirements.  The  graphs  show  that,  as  the  withdrawal  increases  to  a level  above  20  percent  ADF.  the 
percent  of  time  that  withdrawals  cannot  be  made,  either  because  of  natural  flow  limitations  or  passbv 
requirements,  or  both,  will  be  between  60  and  150  days  per  year.  Streams  underlain  by  large  amounts  of 
limestone  are  exceptions  because  they  have  verv'  substantial  base  flows. 
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Figure  2.2.  Impact  of  Selected  Withdrawal  and  Passby  Flow  Combinations,  Ridge  and  Valley 
Freestone,  Wild  Brown  and  Combined  Species 
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Figure  2.3.  Impact  of  Selected  Withdrawal  and  Passby  Flow  Combinations,  Ridge  and  Valley 
Freestone,  Wild  Brook  Trout 
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Figure  2. 4.  Impact  of  Selected  Withdrmval  and  Passby  Flow  Combinations,  Ridge  and  Valley 
Limestone  Group  1,  Wild  Brown  and  Combined  Species 
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Figure  2.5.  Impact  of  Selected  Withdrawal  and  Passby  Flow  Combinations,  Ridge  and  Valley 
Limestone  Group  2,  Wild  Brown  and  Combined  Species 
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Figure  2. 6.  Impact  of  Selected  Withdrawal  and  Passby  Flow  Combinations,  Ridge  and  Valley 
Limestone  Group  1,  Wild  Brook  Trout 
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Figure  2.  7.  Impact  of  Selected  Withdrmval  and  Passby  Flow  Combinations,  Ridge  and  Valiev 
Limestone  Group  2,  Wild  Brook  Trout 
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Figure  2.8.  Impact  of  Selected  Withdrawal  and  Passby  Combinations,  Unglaciated  Plateau 
Segment  Class  1 Streams,  Wild  Brook  Trout 
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Figure  2.9.  Impact  of  Selected  Withdrawal  and  Passby  Combinations,  Unglaciated  Plateau 
Segment  Class  2 Streams,  Wild  Brook  Trout 
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Figure  2. 10.  Impact  of  Selected  Withdrawal  and  Passby  Flow  Combinations,  Unglaciated  Plateau 
Segment  Class  1 Streams,  Wild  Brown  and  Combined  Species 
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Figure  2.11.  Impact  of  Selected  Withdrawal  and  Passby  Flow  Combinations,  Unglaciated  Plateau 
Segment  Class  2 Streams,  Wild  Brcrwn  and  Combined  Species 
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These  impact  graphs  can  be  used  to  develop  statewide  policies  regarding  which  impact  curve  should  be 
used  to  establish  passby  flows.  These  graphs  also  can  be  used  to  determine  impact  of  a proposed 
withdrawal  at  any  site  within  these  stud\'  regions.  These  graphs  also  can  be  used  by  water  purveyors  to 
anaK'ze  stream  intake  alternatives  that  meet  state  fishery  protection  levels  on  cold  water  streams  having 
used  by  water  purveyors  to  anahze  stream  intake  alternatives  that  meet  state  fishery  protection  levels  on 
cold  water  streams  havmg  drainage  areas  less  than  100  square  miles.  The  determination  of  which  impact 
curve(s)  to  use  will  have  to  consider  costs  both  to  the  environment  and  to  withdrawal  users. 

Although  regional  cnteria  have  been  developed,  the  computer  program(s)  can  be  used  to  investigate 
alternatives  or  special  situations  that  have  not  been  considered  in  developing  the  regional  criteria. 
Additional  runs  will  require  hydrology  for  the  study  site(s),  or  for  a stream  where  a withdrawal  is  proposed. 
A regional  hydrology  procedure  has  been  developed  for  use  in  developing  ADF  and  median  monthly  flow 
time  senes  for  any  location  within  these  study  regions  in  Pennsylvania. 

The  regional  hydrology  procedure  uses  data  for  one  or  more  selected  U.S.  Geological  Survey 
stream  gages.  The  hydrologic  regions  are  shown  in  Plate  2.  The  gage(s)  recommended  for  each  hydrologic 
region  are  shown  in  Table  2.1.  The  hydrology'  for  a specific  site  is  generally  derived  by  multiplying  the 
gage  data  by  the  ratio  of  drainage  area  at  the  site  to  drainage  area  at  the  gage.  More  complex  procedures 
are  necessary^  for  sites  with  mixed  or  unusual  geology'.  Adjustments  for  special  conditions  such  as  springs, 
withdrawals,  and  wastewater  treatment  plant  return  flows  are  included  in  the  procedures. 


An  alternative  procedure,  called  the  associated  habitat  duration  impact  analysis  procedure,  is 
designed  to  estimate  the  effect  of  withdrawals  on  the  study  streams  withm  a study  region  or  group.  The 
concept  is  that  impacts  of  a given  level  of  withdrawal  on  a given  study  species  should  be  similar  within  a 
study  region  or  group.  The  procedure  combines  daily  flow  duration  analyses  for  each  season  with  the 
RMWUA  versus  flow  relationship  for  each  study  stream  to  obtam  the  habitat  probability  relationship  for 
each  stream  and  season.  The  habitat  probability  relationships  are  derived  for  unimpacted  and  impacted 
conditions  and  used  to  develop  the  relationship  between  percentage  change  in  habitat  and  flow  for  different 
levels  of  withdrawal.  For  a given  level  of  withdrawal,  the  passby  flow  equals  the  lowest  level  of  flow  for 
which  habitat  reduction  is  equal  to  an  acceptable  level.  Flows,  withdrawals,  and  passby  flows  are 
expressed  as  percentages  of  ADF,  so  that  levels  of  impact  for  different  passby  flows  can  be  compared 
across  streams  within  a study  region. 

These  relationships  can  be  used  to  evaluate  the  effect  of  alternative  withdrawals  and  passby  flows 
on  habitat  over  the  entire  range  of  daily  flows  for  each  season.  The  effect  of  passby  flows  on  water  supply 
availability'  can  be  evaluated  in  a tradeoff  analysis.  This  information  can  be  used  to  determine  allowable 
levels  of  withdrawal,  and  the  required  passby  flows,  to  prevent  unacceptable  impacts  on  the  evaluation 
species.  The  results  can  be  used  to  develop  regional  criteria  for  required  passby  flows.  If  this  procedure 
were  used,  it  could  allow  evaluation  of  the  existing  regional  classes,  which  may  lead  to  modification  of 
those  classes.  The  procedure  is  incomplete  due  to  time  and  cost  constraints  and  the  decision  to  use  the  time 
senes  impact  analysis  procedure  instead. 
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Table  2. 1.  Hydrology  Re}>ions  and  (iaffes 
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Table  2. 1.  Hyilroloi^y  Regions  and  (iai’cs — Continued 
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3.0  STUDY  CONCLUSIONS 


A procedure  has  been  developed  for  determmmg  instream  flow  needs  and  passb\  flows  for  small 
reproducmg  trout  streams  in  Pennsylvania  and  Maiy  land.  The  procedure  is  based  on  available  habitat,  is 
easily  denved  from  hydrologic  records,  and  does  not  require  stream-specific  impact  analysis  studies.  At 
present,  the  procedure  can  be  applied  to  sites  with  drainage  areas  less  than  100  square  miles  in  the  Ridae 
and  Valley,  and  unglaciated  parts  of  the  Appalachian  Plateaus,  physiographic  provinces.  The  procedure 
includes  computer  program(s)  that  estimate  the  impact  on  fisheiy-  habitat  available,  resulting  from  vanous 
combinations  of  withdrawal  and  passby  flow,  for  project  sites  in  those  study  regions.  The  effects  of 
imposmg  passby  flows  on  the  availability  of  water  supply  also  is  estimated.  This  information  can  be  used 
to  evaluate  tradeoffs  between  impacts  on  fishery-  habitat  and  impacts  on  the  water  supply. 

The  computer  program  has  been  used  to  develop  a set  of  graphs  relating  withdrawal,  passby  flow, 
and  impact  on  habitat  for  brook,  browTt,  and  combined  brook  and  brown  trout.  The  impact  of  passby  flows 
on  water  supply  availabiliU'  has  been  superimposed  on  the  habitat  impact  graph  to  facilitate  tradeoff 
analysis  and  development  of  regional  criteria  for  passby  flows.  The  computer  program(s)  also  may  be  used 
to  study  special  situations  not  considered  in  development  of  the  impact  curves.  The  procedures  can  be 
extended  to  the  remaining  parts  of  Pennsylvania,  Maryland,  and  the  Susquehamia  basin  by  collecting  and 
anaK'zing  additional  field  data  for  each  remaining  study  region. 

The  PHABSIM  components  of  IFIM  can  be  applied  to  selected  study  streams  to  develop  the  WUA 
relationships  necessary-  to  estimate  the  impact  of  withdrawals  for  streams  in  a defined  study  region,  and  to 
develop  regional  habitat  impact  curves. 

The  computer  program  developed  as  part  of  this  study  can  be  used  to  determine  the  impacts  of 
withdrawals  for  the  study  sites,  and  the  results  can  be  used  to  develop  regional  relationships  between 
withdrawal,  passby  flow,  and  impact  on  fishery  habitat.  These  relationships  can  be  used  to  develop 
regional  and  statewide  passby  flow  criteria. 

The  original  concept  of  classify-ing  streams  based  on  differences  m key  physical  characteristics  that 
affect  the  availability  of  habitat  at  different  flows  is  safisfactoiy-  for  developing  a regional  procedure  for 
determining  instream  flow  needs. 

The  stream  classification  scheme,  based  on  physiographic  provinces  and  sections,  type  of  geologv-, 
and  stream  segment  number,  appears  to  represent  the  differences  in  the  key  physical  features  that  affect  the 
availability  of  habitat.  However,  the  impact  curves  show  that  there  are  differences  m impact  between 
brook  and  brown  trout,  and  between  brook  trout  and  combined  brook  and  brown  trout.  This  result 
indicates  that  the  trout  species  present  is  an  important  vanable  in  determining  statewide  policy  regarding 
passby  flows. 

The  classification  by  segment  number  is  useful  for  separating  the  impacts  of  withdraw  als  on  small, 
steep  streams  from  those  that  are  larger  and  less  steep.  It  also  is  useful  in  ensunng  that  streams  of  different 
size  are  sampled. 

The  impact  analysis  results  show  differences  in  impacts  between  study  sites  in  different  segment 
classes  in  all  study  regions.  These  differences  are  considered  insignificant  for  the  Ridge  and  Valle> 
Freestone  study  region,  and  impact  curves  for  segment  classes  1,  2.  and  3 were  combined.  Streams  in  the 
Ridge  and  Valley  Limestone  study  region  need  to  be  further  classified  based  on  amount  of  limestone.  The 
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habitat  impact  curves  for  different  segment  classes  behave  erratically,  probably  due  to  site-specific 
differences  in  hydrology,  and  small  sample  size  for  segment  classes  2 through  4.  For  the  Unglaciated 
Plateau  study  region,  the  habitat  impact  curves  are  different  for  sites  in  segment  classes  1 and  2. 


4.0  STUDY  RECOMMENDATIONS 

The  habitat  and  withdrawal  impact  curves  developed  in  this  study  should  be  used  by  the 
participating  agencies  to  develop  regional  or  statewide  procedures  for  determining  withdrawal  limits  and 
passby  flows.  In  particular,  decisions  need  to  be  made  regarding  acceptable  levels  of  impact  on  both  uses. 

This  procedure  also  should  be  extended  to  trout  streams  in  the  Piedmont  Province.  Based  on 
present  knowledge,  it  is  recommended  that  the  province  be  divided  into  the  Piedmont  Upland,  Piedmont 
Lowland  and  GetU  sburg -Newark  Lowland  sections,  and  that  both  limestone  and  freestone  subdivisions  of 
these  sections  be  considered.  Alternatively,  the  entire  province  could  be  classified  as  either  limestone  or 
freestone,  regardless  of  the  physiographic  section. 

The  method  should  be  developed  for  trout  streams  in  the  glaciated  sections  of  the  Appalachian 
Plateaus  Province.  Based  on  present  knowledge,  three  study  regions  are  recommended:  Glaciated  Low 
Plateau  and  Glaciated  Pocono  Plateau  combined;  Glaciated  High  Plateau;  and  Glaciated  Pittsburgh 
Plateau.  Also,  the  study  design  needs  to  consider  the  possibility  that  headwater  streams  formed  on  glacial 
till  are  much  steeper  and  have  different  hydrology  and  habitat  impact  characteristics  than  streams  formed 
on  glacial  fill  materials  in  the  valleys. 

Studies  of  additional  regions  and  t\pes  of  streams  should  include  evaluations  of  the  transferabilitv 
of  HSC  to  these  regions  and  ffpes  of  streams. 

It  has  been  demonstrated  that  regional  relationships  for  fishery  habitat  can  be  developed  for 
Pennsylvania  and  the  Susquehanna  River  Basin  streams.  It  is  appropriate  to  see  if  these  concepts  can  be 
extended  to  larger  cold  water  and  warm  water  streams  and  rivers  in  the  Susquehanna  basin  and 
Pennsylvania.  These  studies  are  needed  because  of  existing  conflicts  between  instream  and  withdrawal 
uses,  and  to  facilitate  evaluation  of  impacts  of  withdrawals  on  those  streams. 

The  applicability  of  results  of  these  studies  to  streams  in  the  Ridge  and  Valley  and  Appalachian 
Plateaus  study  regions  in  Maryland  should  be  considered. 


5.0  AREAS  FOR  ADDITIONAL  RESEARCH 

The  computer  program  should  be  further  refined.  In  particular,  the  hydrology  calculations  that  are 
presently  made  externally  should  be  incorporated  in  the  program.  Also,  a reservoir  operations  model 
should  be  added  to  the  program  to  allow  consideration  of  minimum  releases  from  storage  facilities. 

The  sampling  scheme  utilized  to  select  study  streams  and  segments  generally  provides  satisfactory' 
results.  However,  the  assumptions  used  in  selecting  a sample  of  streams  should  be  investigated  further. 
The  number  of  segment  class  1 study  sites  sampled  appears  to  be  adequate  in  all  study  regions.  The 
number  of  segment  class  2 sites  appears  to  be  adequate  in  the  Ridge  and  Valley  Freestone  and  Unglaciated 
Plateau  study  regions,  but  appears  inadequate  in  the  Ridge  and  Valley  Limestone  study  region.  The 
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number  of  segment  class  3 and  4 sites  appears  to  be  inadequate  in  all  study  regions.  There  ma\-  be  a need 
for  additional  segment  class  3 and  4 study  sites  in  all  study  regions,  and  additional  segment  class  2 sites  m 
the  Ridge  and  Valley  Limestone  region.  Also,  the  relationship  of  the  stream  selection  procedures  to 
vanations  m hydrology  within  a study  region  should  be  evaluated  to  determine  whether  each  hc  drologic 
region  should  be  sampled.  Vanations  in  hydrology  among  segment  classes,  due  to  both  natural  and  man- 
made conditions,  also  should  be  considered. 

Transects  located  near  the  midpoint  of  each  mesohabitat  t\pe  appear  to  provide  satisfactor\- 
sampling  of  spawning  habitat.  In  future  studies,  it  may  be  desirable  to  collect  data  at  a transect  m the 
downstream  part  of  pools  to  include  the  area  with  the  highest  proportion  of  redds. 

The  field  measurement  and  model  calibration  problems  encountered  in  this  stud\-  should  be 
considered  and  minimized  m selecting  streams  for  future  studies. 

The  HSC  developed  in  this  study  are  based  on  the  best  field  data  obtainable  with  the  resources 
available  for  the  study.  However,  these  criteria  could  be  refined  in  future  studies  by;  testing  the  HSC 
developed  in  this  study  against  independent  habitat  usability'  data  for  streams  in  the  same  study  regions; 
developing  separate  HSC  for  each  study  region;  developing  HSC  for  rainbow  trout;  or  collecting  additional 
data  to  allow  evaluation  of  the  effects  of  season,  time  of  day,  or  other  trout  species  present.  Development 
of  habitat  suitability  criteria  for  rainbow  trout  allows  application  of  the  procedures,  including  habitat 
impact  curve  development,  to  that  species. 

The  regional  hydrology  procedures  developed  in  this  study  are  the  best  that  could  be  developed 
within  the  time  and  cost  constraints  of  the  study.  As  expenence  is  gamed  with  the  procedures,  refinements 
may  become  necessary'  or  desirable. 

The  habitat  data  for  the  Mary'land  study  streams  should  be  used  cautiously,  because  of  evidence 
that  some  of  the  streams  are  not  in  dynamic  equilibnum.  The  existing  data  should  be  verified  through  other 
sources,  or  collection  of  additional  data.  Also,  the  effect  of  changes  m bed  and  banks  on  habitat  estimation 
should  be  evaluated. 
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A 


Adult  life  stage 
Average  daily  flow 

ADF 

Associated  habitat 
duration  analysis 

Binary  suitability 
criteria 

Conservation  flow 

Conservation  releases 

Consumptive  use 

Cover 

Daily  flow 
Duration  analysis 


Trout  6 or  more  mches  long. 

The  anthmetic  mean  of  induidual  daily  mean  discharges  during  a period 
of  record. 

Average  daily  flow. 


Development  of  habitat  probability  relationship  b\’  determining  habitat 
corresponding  to  a flow  and  assigning  the  probabilit\'  of  the  flow  to  the 
habitat. 


B 


Habitat  suitabilitv'  critena  that  have  values  only  of  zero  or  unit>-. 


c 

Mandated  flow  expected  to  be  maintained  downstream  from  a water 
storage  facilitv-  or  water  mtake  to  protect  mstream  uses,  including  fisherv 
habitat. 

Releases  made  from  a controlled  water  storage  facility-  to  maintain  some 
amount  of  flow  in  the  stream  downstream  from  the  facilitv'. 

Loss  of  water  from  ground-water  or  surface  water  source,  through  a man- 
made conveyance  system,  by  a process  that  does  not  return  the  water  to 
the  basin. 

Areas  of  shelter  that  provide  resting  places,  visual  isolation,  or  protection 
from  predators  for  aquatic  organisms. 


D 

Average  of  instantaneous  discharges  dunng  a clock  day. 

Categonzation  of  events  (e.g..  flow  rates  or  habitat  available)  to  determine 
the  probabilitv'  of  exceedance  by  arranging  the  values  in  order  of 
magnitude. 


Detailed  analysis  program 


Detailed  analysis  program 

Computer  program  written  in  Microsoft  Excel  format  for  complete 
analysis  of  the  impact  of  any  combination  of  withdrawal  and  passby  flows 
on  the  flow  and  habitat  of  a project  stream;  see  preliminary  analysis 
program. 

E 

Evaluation  species 

Species  used  to  estimate  effects  of  changes  in  flow  on  the  aquatic 
ecosystem. 

Flow  duration  analysis 

F 

Duration  analysis  of  streamflow  data  of  a selected  time  step  (e.g.,  daily  or 
monthly). 

Freestone 

A general  term  for  the  class  of  rocks  that  do  not  contain  significant 
amounts  of  carbonate  minerals.  See  limestone. 

Flow  protection 

Maintenance  of  flows  to  prevent  significant  reductions  in  habitat  for 
aquatic  species,  or  other  instream  uses. 

Fry  life  stage 

Immature  fish  after  emergence  from  gravel,  assumed  herein  to  be  less  than 
2 inches  long. 

Gaging  station 

G 

Point  on  a stream  or  water  body  where  water  surface  elevations  or  flow  are 
systematically  measured. 

Glacial  boundary 

Location  of  the  terminal  moraine  of  the  late  Wisconsin  glacial  advance,  as 
defined  by  Sevon  (1995). 

Habitat 

H 

The  place  where  an  organism  or  population  lives  and  its  surroundings, 
both  living  and  nonliving;  used  herein  to  refer  to  the  physical  aspects  of 
habitat  represented  as  weighted  usable  area. 

Habitat  suitability 
criteria 

Relationship(s)  describmg  usability  of  different  values  of  physical  habitat 
vanable(s)  (depth,  velocity,  substrate/cover)  that  compose  the  physical 
habitat  of  species. 
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HSC 

Habitat  suitability  critena. 

Hydrologic  region 

A portion  of  a study  region  assumed  to  be  h}drologicall\  similar  for 
computing  ADF  and  median  monthly  flows  for  project  streams. 

Habitat  duration 

Duration  analysis  of  habitat  data  of  selected  time  step  (e.g..  daily  or 
monthly). 

Impact 

I 

Absolute  or  percentage  difference  between  the  amount  of  habitat  av  ailable 
without  the  withdrawal  and  the  amount  available  with  the  withdrawal. 

Instream  use 

Any  use  of  water  that  does  not  require  diversion  or  w ithdraw  al  from  the 
natural  watercourse. 

Instream  Flow 

Incremental  Methodology 

A method  to  quantify  the  effects  of  alterations  of  streamflow  on  the 
aquatic  ecosystem. 

IFIM 

Instream  Flow  Incremental  Methodology. 

Inflection  point 

Point  where  the  slope  of  a curve  changes. 

Juvenile  life  stage 

J 

Immature  fish  larger  than  ffy';  assumed  herem  to  be  between  2 and 

6 inches  long. 

Life  stage 

L 

An  arbitraty’  age  classification  of  an  organism  used  in  this  studv’  to 
describe  adult,  juvenile,  ffy-  and  spawning  penods  in  the  life  of  selected 
species. 

Limestone 

A general  term  for  the  class  of  rocks  that  contain  carbonate  minerals 
(calcium  carbonate  or  magnesium  carbonate),  as  shown  by  Pa.  DER 
(1990). 

Limestone  streams 

Streams  draining  areas  underlain  by  carbonate  rocks;  defined  in  this  studv' 
as  streams  having  total  alkalimty  greater  than  70  mg/1,  or  identified  as 
limestone  streams  by  Shaffer  (1991). 
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M 

Median  monthly  flow 

Median  value  of  all  the  daily  flows  during  a particular  month  for  some 
penod-of-record. 

Median  monthly  habitat 

Habitat  available  half  the  tune  during  a particular  month  in  the  record; 
defined  in  this  study  as  habitat  available  at  the  median  monthly  flow. 

Mesohabitat 

Collective  term  for  different  stream  habitat  types  (e.g.,  riffle,  run,  pool). 

Microhabitat 

Small  localized  areas  within  a mesohabitat  type,  typically  described  by  a 
combination  of  depth,  velocity,  substrate,  or  cover. 

No-loss  of  habitat 

N 

No  reduction  of  weighted  usable  area  at  any  flow. 

No-net-loss  of  habitat 

No  reduction  in  weighted  usable  area  at  the  median  monthly  flow. 

No-net-loss  flow 

The  flow  that  results  in  no-net-loss  of  habitat,  computed  as  the  smaller  of 
the  flow  at  the  maximum  renormalized  minimum  weighted  usable  area  and 
the  median  monthly  flow. 

Passby  flow 

P 

The  flow  rate  below  which  a withdrawal  can  not  be  allowed. 

PHABSIM 

Physical  Habitat  Simulation  Program;  a set  of  software  and  methods  used 
to  compute  relationships  between  physical  habitat  and  streamflow. 

Physiographic  province 

Region  with  similar  structural  characteristics  and  a unified  geomorphic 
history,  as  described  by  Fenneman  (1938)  and  delineated  by  Pa.  DER 
(1989)  and  Sevon  (1995). 

Physiographic  section 

A subdivision  of  a physiographic  province,  as  delineated  by  Pa.  DER 
(1989),  Sevon  (1995),  or  Sevon  (in  preparation). 

Pool 

Part  of  a stream  where  velocity  is  reduced,  usually  with  deeper  water  than 
surrounding  areas. 

Preliminary  analysis 
program 

Computer  program  written  in  Microsoft  Excel  format  for  initial  analysis 
of  the  impact  of  combinations  of  withdrawal  and  pre-specified  passby 
flows  on  the  flow  and  habitat  of  a project  stream;  see  detailed  analysis 
program. 
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Protection 

Maintenance  or  protection  of  habitat. 

Recid 

R 

A depression  in  the  streambed  created  by  trout  or  salmon  for  spawning 
purposes. 

Renormalized  minimum 
weighted  usable  area 

The  amount  of  weighted  usable  area  available  for  the  most  limited  life 
stage  at  each  flow,  rescaled  to  a range  of  zero  to  unit\'. 

Reproducing  trout  stream 

Stream  with  naturally  reproducing  trout  population(s). 

Riffle 

Shallow  rapids  in  a stream  where  obstructions  create  waves. 

RiVTVVUA 

Renormalized  minimum  weighted  usable  area. 

Run 

A part  of  a stream  charactenzed  by  rapid  velocit\'  and  few  waves  o\  er  a 
sigmficant  length. 

Season 

s 

Penod  of  time  when  the  same  life  stages  are  present. 

Segment 

A certain  length  of  a study  stream. 

Simulation  flow 

Any  flow  rate  for  which  depth,  velocity'  and  weighted  usable  area  ha\e 
been  computed. 

Spawning  life  stage 

Life  stage  defined  herein  as  including  redd  construction,  lading  and 
incubation  of  eggs,  and  immature  trout  up  to  the  time  of  emergence  from 
the  substrate  in  the  spring  of  the  year. 

Study  region 

A part  of  a physiographic  province  or  section  assumed  to  have 
homogeneous  topographic,  geologic,  hydrologic,  and  habitat 
characteristics. 

Study  site 

A representative  portion  of  a study  segment  selected  for  detailed  data 
collection  and  modeling. 

Study  stream 

A stream  selected  from  lists  of  trout  streams  and  assumed  to  be 
representative  of  other  trout  streams  in  the  same  stud>-  region. 

Substrate 

The  material  on  the  bottom  of  the  stream  channel  such  as  rocks.  gra\  el.  or 
sand. 
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Summer  season 

Months  of  July  through  September,  when  only  adult  and  juvenile  life 
stages  are  present. 

Time  series 

T 

A set  of  values  arranged  in  chronological  order. 

T ransect 

A vertical  cross  section  taken  across  the  stream. 

Univariate  suitability 
criteria 

u 

Habitat  suitability  criteria  that  vary  continuously  over  the  range  from 
zero  to  unity. 

Weighted  usable 

area 

w 

Unit  of  measurement  of  habitat  used  in  Instream  Flow  Incremental 
Methodology;  the  wetted  area  of  a stream  weighted  by  its  suitability  for 
use  by  aquatic  organisms  or  recreational  activity  (units  of  square  feet  per 
thousand  feet  of  stream). 

Wetted  perimeter 

The  length  along  the  bottom  and  sides  of  a stream  channel,  perpendicular 
to  the  flow,  that  is  in  contact  with  the  water  at  a particular  flow  rate. 

Wetted  perimeter 
method 

A method  for  determining  flows  that  maintain  the  availability  of  food 
based  on  the  relationship  of  wetted  perimeter  to  flow. 

WUA 

Weighted  usable  area. 

WWTP 

Wastewater  treatment  plant. 
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Key  to  Study  Sites  Shown  on  Plate  1 


Stream  Name 

Number 

Bear  Run 

1 

Big  Fill  Run.  Seg.  1 

2 

Big  Fill  Run.  Seg.  2 

J 

Big  Run 

4 

Fowler  Hollow.  Seg.  1 

6 

Fowler  Hollow,  Seg. 2 

7 

Green  Creek,  Seg.  1 

9 

Green  Creek.  Seg.  2 

10 

Green  Creek.  Seg.  3 

11 

Homing  Run 

12 

Kansas  Valiev  Run 

13 

Laurel  Run  (Juniata) 

15 

Mile  Run 

16 

Mugser  Run,  Seg.  1 

17 

Mugser  Run,  Seg.  2 

18 

Rapid  Run.  Seg.  1 

19 

Rapid  Run.  Seg.  2 

20 

Rapid  Run.  Seg.  3 

21 

Salem  Creek 

22 

Sand  Spring  Run 

23 

Swift  Run 

24 

Vanscoyoc  Run 

26 

Wapwallopen  Creek.  Seg.  1 

27 

Wapwallopen  Creek.  Seg.  2 

28 

Wapwallopen  Creek,  Seg.  3 

29 

Wapwallopen  Creek.  Seg.  4 

30 

Antes  Creek 

31 

Big  Spring  Creek 

32 

Boiling  Spring  Run 

33 

Bushkill  Creek.  Seg.  1 

34 

Bushkill  Creek,  Seg.  2 

35 

Cedar  Creek  (Lehigh) 

36 

Cedar  Run  (Centre) 

37 

Cedar  Run  (Cumberland) 

38 

Falling  Spring  Run 

39 

Honey  Creek 

40 

Letort  Creek.  Seg.  1 

41 

Letort  Creek,  Seg.  2 

42 

Lick  Creek 

43 

Little  Fishing  Creek 

44 

Long  Hollow  Run 

45 

Monocacy  Creek.  Seg.  1 

46 

Monocacy  Creek,  Seg.  2 

47 

Monocacy  Creek.  Seg.  3 

48 

Nancy  Run 

49 

Penns  Creek,  Seg.  1 

50 

Penns  Creek.  Seg.  2 

51 

Penns  Creek.  Seg.  3 

52 

Potter  Creek 

53 

Stream  Name  | 

Number 

Spnng  Creek  (Berks) 

54 

Spring  Creek.  Seg.  1 

55 

Spring  Creek.  Seg.  2 

56 

Spring  Creek.  Seg.  3 

57 

Spring  Creek.  Seg.  4 

58 

Tnndle  Spring  Run 

59 

Trout  Creek 

60 

Beech  Run 

61 

Benner  Run 

62 

Bloomster  Hollow 

63 

Cherry  Run 

64 

Coke  Oven  Hollow 

65 

Cush  Creek.  Seg.  1 

66 

Cush  Creek.  Seg.  2 

67 

Dunlap  Run 

68 

E.  Br.  Spring  Creek.  Seg.2 

70 

Fall  Creek,  Seg.  1 

71 

Fall  Creek,  Seg.  2 

72 

Findley  Run 

73 

Lower  Two  Mile  Run,  Seg.  1 

74 

Lower  Two  Mile  Run.  Seg.  2 

75 

Lyman  Run 

76 

McClintock  Run 

77 

McEwen  Run 

78 

Meyers  Run 

79 

Mill  Run 

80 

Red  Run 

82 

Seaton  Run 

83 

Strange  Hollow' 

84 

Tannery’  Hollow 

85 

Warner  Brook 

86 

Whites  Creek,  Seg.  1 

88 

Whites  Creek.  Seg.  2 

89 

E.  Br.  Raven  Creek 

90 

Granville  Run 

91 

Laurel  Run  (Huntingdon) 

92 

Baisman  Run 

93 

Basin  Run.  Seg.  1 

94 

Basin  Run.  Seg.  2 

95 

Cooks  Branch 

96 

First  Mine  Branch 

97 

Gillis  Falls.  Seg.  1 

98 

Gillis  Falls.  Seg.  2 

99 

Greene  Branch 

100 

Norris  Run 

101 

Piney  Run 

102 

Third  Mine  Branch 

103 

Timber  Run 

104 

Warren 
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